Electron Spin Relaxation under Drift in GaAs 
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Based on a Monte Carlo method, we investigate the influ- 
ence of transport conditions on the electron spin relaxation in 
GaAs. The decay of initial electron spin polarization is calcu- 
lated as a function of distance under the presence of moderate 
drift fields and/or non-zero injection energies. For relatively 
low fields (a couple of kV/cm), a substantial amount of spin 
polarization is preserved for several microns at 300 K. How- 
ever, it is also found that the spin relaxation rate increases 
rapidly with the drift field, scaling as the square of the elec- 
tron wavevector in the direction of the field. When the elec- 
trons are injected with a high energy, a pronounced decrease 
is observed in the spin relaxation length due to an initial 
increase in the spin precession frequency. Hence, high- field 
or high-energy transport conditions may not be desirable for 
spin-based devices. 

PACS numbers 72.20.Ht, 72.25.Dc, 72.25.Dc, 72.25.Rb 



Recently, there has been an increasing interest in 
the emerging field of spin electronics (spintronics) , 
where the electron's spin degree of freedom is ex- 
ploited. 1-3 Many devices and applications, including 
giant-magnetoresistive structures 4 and magnetic random 
access memories, 5 have been suggested with varying de- 
grees of success. Furthermore, the suggestion of utiliz- 
ing electron spin for quantum information processing is 
a natural extension of spintronics since the electron pro- 
vides an ideal condition for a qubit. 

Of the potential applications, the hybrid devices 6 that 
combine " traditional" semiconductor electronics with the 
utilization of the spin state are currently at the center of 
attention for their increased functionality and the ease 
of integration. Electron spin states in semiconductor 
structures relax (depolarize) by scattering with imper- 
fections or elementary excitations such as other carriers 
and phonons. Therefore, to realize any useful spintronic 
devices, it is essential to understand and have control 
over spin relaxation such that the information is not lost 
before a required operation is completed. So far, most 
of the studies on spin relaxation have been focused on 
electrons with a thermal or near-thermal distribution. 7 
However, electron spins in the spin-based devices can be 
subject to highly non-thermal transport conditions in- 
cluding high drift fields for high-speed transfer of spin 
information. 

In this work, we investigate the influence of transport 
conditions on electron spin relaxation in n-type GaAs. 
A Monte Carlo approach is used to simulate electron 
transport, including the evolution of spin polarization, 8 
to determine spin relaxation lengths and times. A non- 



parabolic energy-band model (with T — L — X valleys) is 
used along with the relevant momentum relaxation mech- 
anisms such as polar optical and acoustic phonon defor- 
mation potential scattering. Of the three main spin relax- 
ation mechanisms [i.e., Bir-Aronov-Pikus, Elliot- Yafet, 
and D'yakonov-Perel' (DP) mechanisms], 7,9 ' 10 we only 
consider the DP processes since this is the dominant 
mechanism in the regime of interest, namely GaAs at 
300K. It is also important to note that the current study 
is limited to the low energy cases, such that electrons are 
found only in the T valley. This is due to the lack of ma- 
terials parameters on spin-orbit spitting in higher energy 
valleys. 

The DP Hamiltonian due to spin-orbital splitting of 
the conduction band may be written as 11 



H - -0-Cl e ff', 



(1) 



and 



^eff - 



■■[k x (ky — k%)x + cyclic perm.] ; (2) 



where a is a dimensionless, material-specific parameter 
which gives the magnitude of the spin-orbit splitting 
a ~ 4rj — rjmo and r\ = A/(E g + A), m is the 

effective mass, k is the electron wave vector, E g is the 
energy separation between the conduction band and va- 
lence band at the T point, and A is the spin-orbit split- 
ting of the valence band. The material parameters for 
GaAs are listed in Table I. 

The quantum mechanical description of the evolution 
of the spin 1/2 is equivalent to the evolution of the clas- 
sical momentum S under an effective magnetic field Q e f f 
with the equation of motion 



-=n eff xs. 



(3) 
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eV 


0.067 


0.0354 


7.0 


1.424 


0.616 


0.341 



TABLE I. Material Parameters for GaAs. m/mo is the 
F valley effective mass ratio, hu>o is the optical phonon en- 
ergy, Ei is the acoustic phonon deformation potential, E g 
is the energy band gap separation at the F point, a np is the 
non-parabolicity factor in the F valley, and A is the spin-orbit 
splitting of the valence band. 
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FIG. 1. Polarization < P > of injected electrons as a func- 
tion of distance. 



Figure 1 depicts the average polarization < P > of 
electrons calculated as a function of distance in GaAs. 
The simulation follows the evolution of 1 x 10 4 spin- 
polarized electrons from an injection plane. Although 
all electrons are initially polarized (< Pq >= 1), polar- 
ization at x = in this figure is not exactly unity due 
to the averaging over a finite mesh size. Clearly, < P > 
decays nearly exponentially with an order of magnitude 
drop (to 0.1) in approximately 4 [im to 7 [ira (5.5 \xm to 
9 fim when the non-parabolicity is not considered) . 

The observed dependence of decay on the applied elec- 
tric field is non-monotonous with the result of 2.5 kV/cm 
placed between the two cases with smaller fields. This 
can be explained by the two competing factors that con- 
tribute to the spin relaxation distance. As the field be- 
comes larger in the linear regime, the electron momen- 
tum and the drift velocity increase in the direction of the 
field. On the other hand, the increased electron momen- 
tum also brings about a stronger effective magnetic field 
as shown in Eq. (2). Consequently, the electron preces- 
sion frequency becomes higher, resulting in faster spin 
relaxation (i.e., shorter spin relaxation time). The inter- 
play of these opposing trends (i.e., higher drift velocity 
and shorter relaxation time) determines the observed be- 
havior of spin relaxation length as a function of electric 
field. 

To examine the dependence of spin relaxation in detail, 
we calculate the spin relaxation rate r" 1 as a function of 
electric field. The value of r s is estimated by approxi- 
mating the spin decay to be exponentially dependent on 
length. The slope of Fig. 1 may be taken to be a charac- 
teristic length Ld, such that 

< P >= Acxp(-L/L D ); (4) 




Electric Field (kV/cm) 

FIG. 2. Spin relaxation rate t^ 1 as a function of electric 
field F. The data points are the results calculated with or 
without the non-parabolicity. The solid line represents the 
< ki > scaling for the non-parabolic band. 



where A is a normalization factor. The spin relaxation 
rate may then be computed as r" 1 — Vdr/Lo, where Vdr 
is the average drift velocity. For comparison, the results 
for the parabolic energy band are shown in Fig. 2 along 
with those with the non-parabolic band. 

As expected, the spin relaxation rate increases rapidly 
with the applied electric field, which is also in accord with 
a recent experimental observation. 12 Although t s for low 
electric fields is nearly the same in both the parabolic 
and non-parabolic band models, the spin relaxation rate 
becomes considerably larger in the non-parabolic case for 
fields higher than approx. 1.5 kV/cm. The values for t s 
at the low fields agree well with those obtained for near 
thermal electrons. 10 

The dependence of spin relaxation on electric field can 
be best understood by considering the functional form 
of the effective magnetic field in Eq. (2). When the ap- 
plied electric field is along the x direction, the interaction 
Hamiltonian reduces approximately to 

H oc (k z z - k y y)k\ . (5) 

Since the average values of k y and k z remain small, it 
is reasonable to assume that r s _1 will scale as < k\ >. 
The solid line in Fig. 2 provides the < k\ > scaling for 
the non-parabolic band, which clearly illustrates good 
agreement over the entire field range under considera- 
tion. This finding indicates that the common assumption 
of electron temperature T to the third power dependence 
[i.e., 1/t<P p cx (/cbT) 3 t p , where t p is the momentum re- 
laxation time] 7 ' 10 substantially overestimates the spin re- 
laxation rates in the drift regime. Considering the nature 
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FIG. 3. Electron polarization < P > vs. distance as a 
function of injection energy and applied electric field. 



of drift in this case, the T 2 scaling may be more appro- 
priate. 

However, the applicability of this scaling rule is also 
limited. As the electron temperature increases, it is clear 
that one needs to take into account the effect of non- 
parabolicity on the spin-orbit splitting (i.e., DP Hamilto- 
nian) itself. By ignoring the higher band contributions, 
Pikus and Titkov obtained the following estimation in 
place of am -3 / 2 in Eq. (2) with the rest of the terms 
unchanged: 7 



can be seen that the electrons quickly lose the extra en- 
ergy and their momentum distribution stablizes within 
1 fim or so, even in the cases of higher energy injec- 
tion (200 meV). Therefore, we have two distinct slopes. 
The first is a more rapid decay in the non-local trans- 
port region, corresponding to a large k and subsequently 
a large precession vector £l e ff. Once thermalization oc- 
curs, the slope becomes smaller and eventually coincides 
with those shown in Fig. 1 with corresponding electric 
fields. For the 200 meV case with an applied field of 1 
kV/cm, the reduction in depolarization length is fairly 
significant, approximately 20 %. 

In summary, we have investigated spin-dependent 
transport in GaAs in the drift regime. Specifically, 
we calculated the spin relaxation time and character- 
istic decay lengths based on a Monte Carlo model. 
The decay of electron spin polarization is found to be 
nearly exponential and the length at which spin relax- 
ation/depolarization occurs was found to be relatively 
large, on the order of 5—10 \im. It is also found that the 
spin relaxation depends strongly on the drift conditions, 
scaling as the square of the electron wavevector in the 
direction of the field. Hence, the commonly assumed T 3 
scaling is not applicable in the drift regime. When the 
electrons are injected with a high-energy, a significant 
reduction is observed in the spin relaxation length. Our 
calculation results are in good agreement with the data 
available in the literature. 

This work was supported by the Office of Naval Re- 
search and the Defense Advanced Research Projects 
Agency. The authors are grateful to Y. Semenov for help- 
ful discussions. 



a(E)m-V 2 (E) = am-V 2 1- 



E{k) 9 - 7r] + 2ry 2 



3-fj 



(6) 



At the constant electric field of 2.5 kV/cm, Eq. (6) re- 
sults in a{E)mT ? '/ 2 {E) that is smaller than am~ 3/2 by 
approximately 10%. The deviation shrinks as the elec- 
tric field decreases, 13 leading to a minor change in the T 2 
scaling of t^ 1 . For the fields beyond the range currently 
under consideration (0.5—2.5 kV/cm), the influence of 
the additional term in Eq. (6) can become crucial. 

In addition to the influence of the drift field, we also 
investigate spin transport and relaxation with varying 
injection energies. As mentioned earlier, one of the fun- 
damental challenges in spintronics has been the electrical 
injection of spin polarized current into a semiconducting 
material. Noting recent suggestions 14 for the utilization 
of tunnel injection, we may anticipate relatively high elec- 
tron energies at the injection point. 

In Fig. 3 we plot the results of simulations where the 
injection energy and applied electric field are changed. It 



1 G. A. Prinz, J. Magn. Magn. Matter 200, 57 (1999). 

2 S. A. Wolf, D. D. Awschalom, R. A. Buhrman, J. M. 
Daughton, S. von Molnar M. L. Roukes, A. Y. Chtchelka- 
nova, and D. M. Treger, Science 294, 1488 (2001). 

3 E. I. Rashba, J. Supercond. 15, 13 (2002). 

4 J. M. Daughton, J. Magn. Magn. Matter 192, 334 (1999). 

5 R. L. Comstock, J. Mater. Sci.-Mater. Electron. 9, 509 
(2002). 

6 See, for example, S. Datta and B. Das, Appl. Phys. Lett. 
56, 665 (1990). 

7 G. E. Pikus and A. N. Titkov, in Opticheskaya Orientat- 
siya (Optical Orientation), edited by F. Meier and B. P. 
Zakharchenya (Nauka, Leningrad, 1989). 

8 A. A. Kiselev and K. W. Kim, Phys. Rev. B 61, 13115 
(2000). 

9 J. Fabian and S. Das Sarma, J. Vac. Sci. Technol. B 17, 
1708 (1999). 

10 P. H. Song and K. W. Kim, Phys. Rev. B 66, 035207 (2002). 

11 M. I. D'yakonov and V. I. Perel, Fiz. Tverd. Tela 
(Leningrad) 13, 3581 (1971) [Sov. Phys. Solid State 13, 



3 



3023 (1972)]. 

H. Sanada, I. Arata, Y. Ohno, K. Ohtani, Z. Chen, K. 
Kayanuma, Y. Oka, F. Matsukura, and H. Ohno, to be 
published in Proc. of the 26th Intl. Conf. on the Physics of 
Semiconductor. 

The effect of non-parabolic energy band in Eq. (6) may be 
non-zero even at very small electric fields due to the finite 
average electron energy. For the average electron energy as 
a function of applied electric field, see K. Tomizawa, Nu- 
merical Simulation of Submicron Semiconductor Devices 
(Artec House, Norwood, MA, 1993), p. 102. 
E. I. Rashba, Phys. Rev. B 62, 116267 (2000). 



4 



